Astronomy & Astrophysics manuscript no. bastian © ESO 2013 

March 1,2013 



Spectroscopic Constraints on tlie Form of the Stellar Cluster Mass 

Function 

N. Bastian^ I.S. Konstantopoulos^ , G. Trancho^-^, D.R. Weisz^, S.S. Larsen^, M. Fouesneau^, C. B. Kaschinski', and 

M. Gieles** 

' Excellence Cluster Universe, Boltzmannstr. 2, 85748 Garching, Germany 
e-mail: bastian@usm.lmu.de 
- - - ~ Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802, USA 

Cn I 3 Gemini Observatory, Casilla 603, La Serena, Chile 

y^ , * Giant Magellan Telescope Organization, 25 IS. Lake ave., Pasadena 91001,USA 

O . 5 Department of Astronomy, Box 351580, University of Washington, Seattle, WA 98195, USA 

CN ' '^ Department of Astrophysics, IMAPP, Radboud University Nijmegen, P.O. Box 9010, 6500 GL Nijmegen, The Netherlands 

$_( I ' Institut fiir Astronomic und Astrophysik der Universitat Miinchen, Scheinerstrafie 1, 81679 Munchen, Germany 

Cy , ** Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 OHA, UK 



(N 



Received Feb. ZZ, 2012; accepted YYY ZZ, 2012 

ABSTRACT 



This contribution addresses the question of whether the initial cluster mass function (ICMF) has a fundamental limit (or truncation) 

^^ ' at high masses. The shape of the ICMF at high masses can be studied using the most massive young (<10 Myr) clusters, however this 

has proven difficult due to low-number statistics. In this contribution we use an alternative method based on the luminosities of the 

brightest clusters, combined with their ages. The advantages are that more clusters can be used and that the ICMF leaves a distinct 

(-H . pattern on the global relation between the cluster luminosity and median age within a population. If a truncation is present, a generic 

Oj' prediction (nearly independent of the cluster disruption law adopted) is that the median age of bright clusters should be younger than 

I I that of fainter clusters. In the case of an non-truncated ICMF, the median age should be independent of cluster luminosity. Here, we 

O ■ present optical spectroscopy of twelve young stellar clusters in the face-on spiral galaxy NGC 2997. The spectra are used to estimate 

the age of each cluster, and the brightness of the clusters is taken from the literature. The observations are compared with the model 

C/3 ' expectations of Larsen (2009) for various ICMF forms and both mass dependent and mass independent cluster disruption. While 

Cu ' there exists some degeneracy between the truncation mass and the amount of mass independent disruption, the observations favour 

a truncated ICMF. For low or modest amounts of mass independent disruption, a truncation mass of 5 - 6 x IO'Mq is estimated, 

consistent with previous determinations. Additionally, we investigate possible truncations in the ICMF in the spiral galaxy M83, 

K^ the interacting Antennae galaxies, and the collection of spiral and dwarf galaxies present in Larsen (2009) based on photometric 

. J^ ' catalogues taken from the literature, and find that all catalogues are consistent with having a truncation in the cluster mass functions. 

iL •« . However for the case of the Antennae, we find a truncation mass of a few x 10* Mq, suggesting a dependence on the environment, as 

1^ . has been previously suggested. 

lO ■ Key words. Galaxies: star clusters: general 

cn 

o 

OQ 1 . Introduction well described by a power-law with index 2, at least on the 

'""' , low mass end (e.g., de Grijs et al. 2003). 

>■ Since their discovery in starburst and merging galaxies (e.g., .^ • ■ • • ■•■••„ • ■ 

c u • inoT II 1. .1 inm\ ° '^ ° • ^ i Dcspite these intensive studies, there is still ongoing de- 

• ^H Schweizer 1987, Holtzman et al. 1992), young massive stel- , . \ ,. r ,\ ■ ■ r ? 

k> 1 1 . /v^. xr^ 1 1 u 1 1 . bate in the literature on some or the basic properties or cluster 

>v . lar clusters (YMCs, a.k.a. young globular clusters or super ,. ,, -■• ^, .i-, 

^^ ^ 1 » X u ^. .u r r 1 .u r . J- / r populations. Is the mass runction or clusters truncated at high 

H ' star clusters) have been the locus ol a plethora ot studies (c.l. , ^. , , ^^^^ , „ ■ -,r,r,o t -.r,r>>. 

. =?.■ Portegies Zwart, McMillan, & Gieles 2010). These studies have ^f "'o^ni.-'^' Z 5n?of' '.'t" ' V"""'™ ' 

c , • J- J 1 1 . / c vu o /-■ 11 u ^rvn1 Gicles 2009; Bastian et al. 2012) or docs the powcf-law continuc 

locused on individual clusters (e.g.. Smith & Gallagher 2001; ■ ■ • ■ , / ^, ■ , -.^,^ ,,„ • 

T . 1 ^rvrvo\ 11 \c 11 i . i ? / to the highcst mass clustcrs (e.g., Chandar ct al. 2010; Whitmore 

Larsen et al. 2008) as well as lull cluster populations (e.g., , ^„f„.nii i i i i- i. i i- 

T^ , , , /, r,f^[^n A J .1 onin\- If et al. 2010) .'How long do clusters live, and how does disruption 

Konstantopoulos et al. 2009; Adamo et al. 2010) in a variety ol ^ , , , • r, , , .• • . . 

, 4. c u 4 ■ »• c^r-\«r^ 1 »• anect the cluster population .' Is cluster disruption dependent on 

galaxy types. Some characteristics ol YMC populations appear , ^ ■ , ■ ■ ■ n ^ ■ i i- 

f , , f , ^ ^, . r, ^, , ^ , -^u . ■ . , the mass ol the cluster and Its environment .' Or IS cluster disrup- 

to be related to the properties ol the host galaxy, with high star- ... , -■ • , • ■ • ■ , • , 

c ,. , ,^.^,. .u.\j- tion independent ol environment, being dominated by internal 

lormation rate galaxies (i.e., starbursts) producing more mas- ^ , _ . ^_,, ^ ° • sr, ,, 

, ^ , ■, ninnn, r> .• ^^v^vo^ j.u processes (sce Bastian 201 1 tor a reccnt review) .' Here, we lo- 

sive clusters (e.g., Larsen 2002, Bastian 2008), and there are ^ , j, ^ , . . . , j. . ^ , 

.... ,u\^u 1 . r .• £c • /• .u r cus on the lorm ol the initial mass lunction ol clusters, paying 

indications that the cluster lormation emciency (i.e., the Irac- . , . ■ ■ • • ■ ■ , • ■ 

t ^ ^u ^ c • 1 . \ • vu • • . particular attention to whether it is truncated at the high mass 

tion ol stars that lorm in clusters) increases with increasing star- ^ , ° 

end 
formation rate surface density of the host galaxy (e.g., Meurer 

et al. 1995; Zepf et al. 1999; Larsen 2004; Goddard et al. 2010; Whether the initial mass function of stellar clusters (ICMF) 

Adamo et al. 201 1, Silva- Villa & Larsen 201 1). Additionally, the has a truncation at the high mass end is important for a number of 

form of the cluster luminosity or mass distributions appears to be reasons. In particular, if there is a truncation, and the truncation 
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varies as a function of environment (e.g., Gieles et al. 2006b; 
Bastian et al. 2012) this may be reflecting important physics in 
the star/cluster formation process, or about the distribution of 
the ISM within galaxies (i.e, size scale of GMCs). Additionally, 
a truncation of the ICMF would have important implications for 
the globular cluster specific frequency problem (e.g., Kruijssen 
& Cooper 2012), with the number of clusters surviving after a 
Hubble time being strongly dependent on whether a truncation 
occurs at ~ 10^ Mq or ~ 10* Mq. Finally, knowing the form 
of the ICMF allows for tighter constraints to be placed on the 
role of disruption in shaping cluster populations (e.g, Gieles & 
Bastian 2008; Chandar et al. 2010; Bastian et al. 2012). 

The form of the ICMF and the mechanism of cluster 
disruption leaves a characteristic imprint on the relation be- 
tween the age and absolute magnitude of ensembles of clusters 
(Larsen 2009). For example, if cluster disruption is strong (i.e., 
the timescale of disruption is short) and independent of environ- 
ment (a mass-independent dissolution time scale, e.g., all clus- 
ters lose 90% of their mass every decade in age - Fall et al. 2009), 
then few old (> 50 Myr) bright clusters are expected to exist 
within the population. Additionally, if the cluster mass function 
is truncated at the upper end (i.e., the mass function is described 
by a Schechter function), this also limits the likelihood of finding 
a bright, old cluster 

Larsen (2009) derived the expected relations between the 
median age vs. absolute magnitude of a cluster population, as- 
suming a range of ICMFs and disruption laws. He used pub- 
lished catalogues of the brightest and 5th brightest clusters in a 
sample of spiral and dwarf galaxies, and found evidence that the 
mass function was indeed truncated at a value of ~ 2 x 10^ Mq, 
and that strong mass independent disruption (90% per decade in 
age) scenarios were not favoured. 

In the present work, we apply the Larsen (2009) analysis 
to a single galaxy, NGC 2997, and use Gemini spectroscopy 
to derive the age of the individual clusters, eliminating the 
age/extinction/metallicity degeneracies known to hamper pho- 
tometric studies. 



2. Observations 

NGC 2997 is a face on spiral galaxy, located at a distance of 
~ 9.5 Mpc (distance modulus of 29.9 - de Vaucouleurs 1979). 
It hosts a relatively rich young cluster population (Larsen & 
Richtler 1999), including one of the brightest YMCs known in a 
spiral galaxy (My ~ -12.9; Larsen 2009). The face on orienta- 
tion of NGC 2997, as well as its relatively large distance, allows 
us to sample the full galaxy with a single pointing of GMOS 
(discussed in more detail in § 12.21 ). mitigating any potential se- 
lection biases of focussing on the central regions of the galaxy, 
giving a more global picture. 

2.1. NGC 2997 cluster sample 

The sample of clusters was taken from the catalogue of 
Larsen (2004), which is based on ground-based U, B, V, Ho- 
and I-band imaging. We refer to Larsen& Richtler (1999) and 
Larsen (1999) for details of the cluster selection process and 
derivation of the completeness limits, and only summarise here. 
Cluster candidates were selected through DAOFIND routine in 
IRAF, and a colour cut of B - V < 0.45 was applied to remove 
foreground stars. Additionally, magnitude cuts were applied 
in order to remove individual bright stars within NGC 2997, 
namely My < -9.5 for blue sources (U - B < -0.4) and My 



< -8 .5 for red sources (U - B > -0.4). Artificial star and cluster 
tests resulted in approximately 80% completeness at My= -9.5. 
All candidates were visually inspected and large diffuse sources 
(likely background galaxies or extended high-surface brightness 
regions within spiral arms) were removed. 

Additionally, sources that had significant Ha emission as- 
sociated with them were removed, which has two benefits. The 
first is that it is often difficult to distinguish a stellar cluster from 
an unbound association at young ages (e.g., Larsen 2004; Gieles 
& Portegies Zwart 2011; Bastian et al. 2011). However, once 
a stellar group has had time to dynamically evolve, it is eas- 
ier to distinguish between clusters and associations (Gieles & 
Portegies Zwart 2011), so restricting to candidates older than 
5-7 Myr helps in removing many spurious sources (i.e., un- 
bound associations). Secondly, it has been suggested that the 
rapid removal of gas from a young cluster can disrupt a large 
fraction of clusters or can lead to significant stellar mass loss 
within clusters (e.g., Goodwin & Bastian 2006; Goodwin 2009). 
However, within a few crossing times, the clusters should sta- 
bilise (e.g., Portegies Zwart et al. 2010), so limiting our sample 
to older clusters lessens the impact that infant mortality or infant 
weight loss may have on the population. 

2.2. Spectroscopy 

The brightest 30 candidates were chosen, and a multi-slit spec- 
troscopic mask was made that maximised the number of clusters 
that could be observed in a single pointing. When the dispersion 
axis of two potential targets overlapped, the brighter cluster was 
chosen. This resulted in 12 clusters that could be observed in a 
single pointing. The absolute V-band magnitudes of these clus- 
ters were taken from the Larsen catalogue and the basic proper- 
ties of the clusters are given in Table [1] 

We obtained spectra of 12 sources in NGC 2997 using the 
Multi-Object Spectroscopy (MOS) mode of the Gemini Multi- 
Object Spectrograph (GMOS) on Gemini South. The data were 
obtained as part of Queue program GS-2009B-Q-29 (PI N. 
Bastian). We used the B600 grating and a slit width of 0.75 
arcsec, resulting in an instrumental resolution of 110 kms"' at 
5000 A. The spectroscopic observations were obtained as four 
individual exposures with an exposure time of 1800 s each (for 
a total on-target integration time of 7200 s). No atmospheric 
dispersion corrector (ADC) was available on Gemini South at 
the time of the observations, causing wavelength-dependent slit 
losses. We corrected each exposure for this effect using the 
method of Filippenko (1982). 

The basic reductions of the data were done using a combi- 
nation of the Gemini IRAF package and custom reduction tech- 
niques, as described in Appendix A in Trancho et al. (2007). 

The age of each cluster was derived using the method pre- 
sented in Trancho et al. (2007). In brief, a template is built for 
each cluster (avoiding lines affected by nebular emission) using 
simple stellar population (SSP) models (Gonzalez Delgado et 
al. 2005) and the PPxF method (Cappellari & Emsellem 2004). 
Spectral indices are measured from the template spectra ( H+He, 
K, m, HyA, M^/75177,Fe5270,Fe5335; see Trancho et al. 2007 
for the full definitions), along with error estimates, which are 
then compared to SSP models (measured in the same way) in 
a least x^ sense in order to choose the best fitting model. The 
uncertainty in age is derived by running 2000 monte carlo sim- 
ulations for each cluster, adding random noise to the measured 
indices, sampled from the estimated errors. The age uncertainty 
is then taken as the standard deviation of a log-normal fit to 
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Fig. 1. The observed spectra of the 12 observed clusters in 
NGC 2997 presented in this work. The ID (following L04) and 
derived spectroscopic (logarithmic) age (in yr) are given for each 
target. 



the resulting age distribution. This method has been shown to 
be largely independent of the signal-to-noise ratio of the data 
(Trancho et al. 2007). 

The observed spectra are shown in Fig. [T] along with the es- 
timated age of each cluster. Additionally, the cluster properties 
are given in Table 1. Only the 12 clusters with spectroscopic age 
determinations are used for determining the cluster mass func- 
tion. 



3. Comparison with IVIodels 

The underlying initial cluster mass function (ICMF) and type of 
disruption strongly affect the relation between the observed age 
and absolute magnitude of a cluster population. Larsen (2009) 
presented an analytic representation of the cfiN jidtdL) distribu- 
tion, i.e., the number of clusters with an age between t and t + dt 
and luminosity between L and L + dL, that considers functional 
forms for the ICMF and the mass evolution of the stellar pop- 
ulation due to disruption. From this it is possible to derive the 
median age (as well as the mean and the standard deviation) as 
a function of My. He considered various ICMF types (truncated 
or untruncated) and different disruption laws (mass/environment 
dependent or independent). 

If disruption does not influence the population, and 
the mass function has an exponential cutoff (i.e., is de- 
scribed by a Schechter (1976) type function - dN/dMi ~ 
Mr^exp(-M,7M*)) then the brightest clusters in a population 
are preferentially young (see also Bastian 2008). The exact re- 
lation between median age and My depends on the truncation 
mass, M*, with lower M* values resulting in the brightest clus- 



ters being younger If cluster disruption is included, in particular 
if cluster disruption is dependent on the cluster mass, as expected 
from N-body and analytic considerations (see e.g., Portegies 
Zwart, McMillan, & Gieles 2010) then the older clusters and 
lower mass clusters are most aff'ected. This results in a turn- 
around in the relation between age and My with fainter clusters 
also being preferentially young. Alternatively, if the ICMF does 
not have a truncation and there is no mass dependent disruption, 
then the median age of a cluster population should be indepen- 
dent of My (the vertical dash-dotted line in the first panel of 
Fig. |2]l. Including strong mass independent disruption (where a 
large and constant fraction of the population disappears every 
decade of ago also results in the median age being indepen- 
dent of My. 

Three characteristic cases are shown in Fig. |2] (after 
Larsen 2009). In each of the panels six models are shown. In 
all models a constant star/cluster formation rate is assumed. The 
dashed lines represent models that adopt a Schechter ICMF with 
M* values of 1 x 10^ 1 x 10^ 2 x 10^1 x 10^ and 1 x IO^Mq, 
from left to right, respectively. In all cases we adopted an index 
of -2 for the ICMF at the low mass end. The dash-dotted lines 
represent the results of models that have a pure power-law ICMF 
(untruncated). The left panel shows models that have no disrup- 
tion, while the middle panel shows models with mass dependent 
disruption (adopting the formalism of Lamers et al. (2005), with 
f4 = 100 Myr), and the right panel adopts no mass dependent dis- 
ruption but instead strong mass independent disruption (a long 
duration infant mortality rate of 90% - see Whitmore et al. 2007). 
The models explicitly include the effects of stochasticity in the 
cluster population and do not require a fully sampled mass func- 
tion at young ages. 

In all cases, the models represent the median age of the clus- 
ters for a given magnitude. Hence if the model is a good repre- 
sentation of the data, half of the clusters should lie to the left and 
half to the right of the curve, at all magnitudes. From Fig. |2]it is 
clear that certain models do not represent the observed data well. 
Low or high truncation values (i.e., ~ lO'* or ~ lO^M©) are not 
favoured if there is little or no mass dependent disruption, and 
low truncation values are also not favoured in the case of strong 
mass independent disruption. 

The data presented in Fig. |2] show on overall tendency for 
younger clusters to be brighter than older ones (Spearman's 
rank correlation coefficient of 0.58, meaning that there is ~ 5% 
chance of log(age) and My to show such a relation by chance). 
This is a generic prediction for any cluster population that dis- 
plays a truncation in the upper end of the mass function (e.g., 
Gieles et al. 2006a, Bastian 2008, Larsen 2009). If no cluster dis- 
ruption is assumed, then the data points are most consistent with 
a truncation mass of a few times lO^M©. Additionally, while the 
impact of mass dependent disruption is largely negligible for the 
observed sample (due to the fact that the brightest clusters tend 
to be young and massive), strong mass independent disruption 
models are not favoured by the data. 



' There is a subtle difference between all clusters losing a fraction 
of their mass every decade in age, and a constant fraction of clusters 
are completely disrupted. If the mass function is a pure untruncated 
power-law (with an index of -2), then there is no difference in terms 
of the method used here. If there is any structure in the ICMF (such as 
a truncation), in the case of all clusters losing a fraction of their mass, 
then this structure will shift towards lower masses with time. 



N. Bastian et al.: Constraining the cluster initial mass function 



No Disruption 



t,= 100 Myr, Fm,„=0 



= 0.9 



518 
51 



■ I 

■;. \^ 679 
!t^jy693 



I \ 
1 \ 

I \ 



g M^ ^^ 4.0 5.d\5.3\6.0\7.q| 
■ ■■*'■■■■*■■■■**■■■■**■■*■■"* 



.1_I_L 



■ < ' 

■ > " 



1 ^ 

I > V I\ 



I I 

> ■ ; ,' 

1 ■ I ,/ 

1 1 # ,/ 

1 11,/ 

1 tit/' 

■I I I J ■»'■ I ■ 



■ \ ^ T ' ; ■' 



,"T' 

'Mil I 



■II, 
Ml, 
'■11 



1111 
I 1 I 



,1— ?I— I 



■ I ■ ■ .1. li.t. .1 ■ II 



7.0 7.5 8.0 8.5 

log Age/yr 



7.0 7.5 8.0 8.5 

log (Age/yr) 



7.5 B.O 

log (Age/yr) 



Fig. 2. Age vs. My for the twelve clusters in the present study. In each panel, we also show the expectations of the median My, 
for a given age, following the results of Larsen (2009). The model disruption parameters are given at the top of each panel, where 
4 is the disruption timescale of a 10"* M© cluster in the mass dependent cluster disruption scenario, and Fmid is the fraction of 
clusters destroyed every decade in age. The dashed lines represent models with different values of the Schechter truncation mass, 
and the dashed-dot lines in each panel represents the case without truncation (i.e., a pure power-law). In the middle and right panels, 
the original data are shown in grey, while the binned data (the mean in flux and age, with the error bars representing the standard 
deviation) are shown as solid (blue) triangles. Two diff'erent bin sizes are shown. 
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09 45 51.34 


-31 10 56.3 
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Table 1. The ID, location, ranking (in terms of V-band absolute brightness in the population), spectroscopicafly determined age and 
error, and absolute V-band magnitude of the clusters in NGC 2997 presented in the current work. 



3.1. Statistical measures 

In order to place the above conclusions on a more statistically 
robust footing, we have also carried out maximum likelihood 
tests, comparing the observations to a wide range of models. For 
the models, we use the Larsen (2009) analytic prescriptions, and 
vary M* from 1 x 10"* to IxIO^Mq, 4 from 100 to 10^ Myr, and 
Fmid from to 0.9. 

To assess the likelihood of our data for a given model, we 
assume that the uncertainties on both age and My follow a 
Gaussian distribution. From this assumption, we can then write 
down the likelihood function to be: 



X.OZ e' 



where 



(1) 



X 



!=1 



(log{age\ots - log(age)i^„,od)^ (Mvj,obs - Mv,i,modY 



(2) 



\og{age)ij,t. 



Mvj.ob, 



and A^ is the number of observed data points. Because we are 
only interested in relative likelihoods, we do not compute the 
normalisation term for the likelihood function in Equation[T] For 
this analysis we limit the comparison between observations and 
models to the range of magnitudes of the observations. 

Using this framework, we sampled the likelihood function 
for 6300 models described in Larsen (2009). The ensemble of 
models that maximise the likelihood function are deemed to rep- 
resent a set of the most likely model parameters (e.g., truncation 
mass, 4, etc). 



3.2. NGC 2997 results 

As expected from visual inspection of Fig. |2] the results of 
the maximum likelihood fitting were largely independent of the 
mass dependent disruption timescale, 4. However, for all values 
of 4, a truncation in the mass function was favoured by the fitting 
procedure. In order tojudge the importance of mass independent 
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disruption (i.e., long duration infant mortality) we set the mass 
dependent disruption timescale to 10^ Myr (i.e., essentially in- 
finite). Doing this, we are left with only two free parameters in 
the model, M* and Fmid- 

The results from this two dimensional maximum likelihood 
analysis are shown in Fig. |3] The darker the shading the higher 
likelihood (also shown in coarser steps as dashed contour lines). 
A clear maximum in this space is seen, a largely horizontal "val- 
ley" with a truncation at 5 x lO^M© that increases slightly as 
a function of Fmid- The solid (white) circles denote the models 
with the maximum likelihood (i.e., with the same likelihood). 
While eleven of the models have M* between 5 - 6 x IO^Mq and 
Fmid:^ 0.2, two of the models had significantly higher Fmid val- 
ues (0.5 and 0.7) and somewhat larger M* values (l-2xlO^M0). 
These models are shown in Fig. |4] 

Hence, we find that a truncation in the cluster mass func- 
tion is favoured independent of the disruption law adopted. If 
disruption is dependent on mass and the surrounding environ- 
ment of a cluster, then a truncation of ~ 5 x lO^M© is strongly 
preferred. A similar conclusion is reached, that a truncation of 
~ 5 - lOx 1O''M0 is required, if the fraction of clusters disrupted 
every decade in, Fmid, is less or equal to 0.5. For higher val- 
ues of Fmid a truncation is still preferred, although at higher M* 
values. 

The largest uncertainty in the value of M* is due to the dis- 
ruption law adopted. If we adopt a mass and environmentally 
dependent disruption scenario, as found for M83 (Bastian et 
al. 2011, 2012) and dwarf galaxies (Cook et al. 2012), then the 
uncertainties are drastically reduced, and instead are due to the 
limitations of age dating clusters, resulting in an uncertainty in 
M* within a factor of two. 



4. Other cluster catalogues 

Additionally, we have run our maximum likelihood method over 
the data presented in Larsen (2009). This catalogue contains 
the brightest and 5* brightest cluster in a photometrically deter- 
mined sample of cluster populations in spiral and dwarf galaxies. 
The fits result in a very similar picture, with M* slightly lower, 
2 X IO^Mq. 

Finally, we can look at the photometrically determined prop- 
erties of cluster populations, taken from the literature. For this, 
we use the 50 brightest clusters in the Antennea galaxies as 
determined by Whitmore et al. (2010) and the sample of clus- 
ters in M83 taken from Bastian et al. (2011). For each dataset, 
we removed sources with ages less than 6 Myr (see § 12.11 for 
a discussion on this point) and then looked at the median age 
and brightness of the population in bins of 4 and 10 clusters 
in the Antennae and M83, respectively. The data are shown in 
Fig. |5] Additionally, for the M83 results, we have overplotted 
the expected model as derived from the full population study of 
Bastian et al. (2012). These are models with M*= 1.6 x IO^Mq 
and 4- 130 Myr for the inner field and M*= 5 x IO'^Mq and 
4- 600 Myr for the outer field. 

Both galaxies are consistent with having a truncated mass 
function, with the Antennae having an M* value of a few 
xlO^Mo, while M* is ~ lO^Moin M83. The value of ~ IO^'Mq 
found for the Antennae galaxies is similar to that suggested by 
Zhang & Fall (1999) and Gieles et al. (2006a) based on the 
observed luminosity function of the cluster population. As ex- 
pected, the cluster mass function extends to higher mass than 
seen in spiral or dwarf galaxies, however it still appears to be 
truncated. 
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Fig. 3. Results from the two dimensional statistical analysis (i.e., 
adopting no mass dependent disruption, and determining the 
likelihood for the models only varying M* and Fmid)- The 
shaded region/contour lines denote the logarithm of the likeli- 
hood of the model for the combination of M* and Fmid, and 
the filled (white) points denote the maximum likelihood models. 
Top: The analysis carried out on the unbinned dataset. Bottom: 
The same analysis, but now carried out on the binned dataset 
(four bins, each with the mean age/flux of three clusters). The 
white points from the upper panel denote models that are shown 
in Fig. H 

5. Conclusions 

We confirm the findings of Larsen (2009) that the cluster mass 
function in spiral galaxies is consistent with a Schechter mass 
function with a truncation near 2 - 5 x IO^Mq. Additionally, 
strong mass independent disruption (> 80 - 90% per decade in 
age) is not supported by the current observations. Comparing the 
age-Mv method to that derived with full population studies in 
M83 leads to largely consistent results (i.e., that a truncation in 
the ICMF is needed and that strong mass independent disruption 
is disfavoured). Applying this method to the cluster population 
of the Antennae galaxies shows that it is also consistent with 
a truncated ICMF, however the truncation mass is significantly 
higher, a few xIO^Mq. 
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Fig. 5. Age vs. My for clusters in the Antennae (left panel - Whitmore et al. 2010), M83 inner region (middle panel - Bastian et 
al. 2011) and M83 outer region (right panel - Bastian et al. 2011). Only clusters with ages greater than 6 Myr were used. In all 
panels the brightest cluster in the sample is shown as a (green) filled circle. For the Antennae, the triangles represent the median 
of the age and My of clusters in the rest of the catalogue (in bins of four clusters). For M83, each bin contains 10 clusters. In the 
middle/right panels, the model corresponding to the best fit results of Bastian et al. (2012) are shown as solid (red) lines. They 
have been corrected the a Salpeter stellar IMF adopted by the Larsen (2009) models in M*. For the Antennae data, the error bars 
represent the standard deviation for each bin. For the M83 data, the error bars represent the region where 60% of the data in that bin 
lie. We note that the models also show such scatter, i.e. only the median of the models is shown. 
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Fig. 4. The observed age vs. My relation along with the maxi- 
mum likelihood models from § 13.21 The solid (blue curves) rep- 
resent models that have Fmid=0, and from left to right have 4 
values ranging from 100 Myr (left) to 10 Gyr (right). All these 
models have M* values between 5 and 6xlO^M0. The dashed 
(red) lines have similar M* values, but no mass dependent dis- 
ruption and a small amount of mass independent disruption. The 
dash-dotted (green) lines also represent models with a truncation 
mass (> IO^Mq) but have higher mass independent disruption 
rates. 
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